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We optimize the signal-to-noise ratio of a Mach-Zehnder atom interferometer with Gaussian 
squeezed input states, in the presence interactions. For weak interactions, our results coincide 
with Phys. Rev. Lett. 100, 250406 (2008), with optimal initial number-variance a oc 7V 1//3 and 
optimal signal-to-noise ratio s a oc TV 2 / 3 for total atom number TV. As the interaction strength 
u increases past unity, phase-diffusion becomes dominant, leading to a transition in the opt imal 
squeezing from initial number-squeezing to initial phase-squeezing with a oc V uN and s oc yj N/u 
shot-noise scaling. The initial phase-squeezing translates into hold-time number-squeezing, which is 
less sensitive to interactions than coherent states and improves s by a factor of y/u. 
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I. INTRODUCTION 

Heiscnbcrg limited atom interferometers offer the pos- 
sibility for compact, inexpensive measurement tools 
which will eventually operate at an unprecedented level of 
precision. Bosc-Einstcin condensates (BECs) in double- 
well potentials are well suited as a platform for such 
devices, as demonstrated by series of recent matter- 
wave interference experiments fH-fi~3| . These experiments 
demonstrate that the double-well BEC system has the 
necessary phase-coherence, and a capacity for fine-tuning 
of the tunneling and interaction parameters, necessary to 
operate an atom interferometer at the maximum sensi- 
tivity allowed by quantum mechanics. 

A doube-well BEC interferometer is typically based 
on the Mach-Zhender Interferometer (MZI) paradigm, in 
which a bimodal input state is mixed by a 50/50 beam- 
splitter, then held for a fixed duration while the two 
modes acquire a relative phase differential 6 (via an ex- 
ternal field), and then mixed again by a second 50/50 
beamsplitter. A measurement of the particle number- 
difference at the output then acts as an estimator for the 
accumulated phase differential. Aside from the prepara- 
tion of the initial state, atom-atom interactions are typ- 
ically neglected in theoretical treatments of the MZI. 

If the input state is a two-mode coherent state (i.e. 
each particle is in the same single-particle orbital), the 
phase-estimation uncertainty, AO is governed by the 
Standard Quantum Limit (SQL), often referred to as the 
"shot-noise-limit", for which AO = l/\/N, with N be- 
ing the total number of particles used. In pioneering 
early work (l3 - [T8j . it became understood that quantum 
mechanics ultimately limits the phase-estimation uncer- 
tainty at the so-called Heisenberg limit Ad = 1/N, which 
is a factor V^V below shot-noise. In a typical realization, 
the best one can do is demonstrate "Heiscnbcrg-Limit 
scaling", which means that AO = q/N, where q is an N- 
independent constant. To reach the Heisenberg limit in 
a bimodal MZI, one must prepare a strongly number- 
sqeezed input state. While the maximally squeezed 
Twin-Fock state (TFS), having exactly N/2 atoms in 



each mode, can achieve Heisenberg scaling at = 0, this 
requires two or three measurements fl9| . while for ^ 0, 
the TFS actually performs worse than shot-noise, and 
should thus be avoided. As an alternative to the TFS, the 
Gaussian Squeezed State (GSS) with optimized squeez- 
ing exhibits a single-measurement phase-estimation un- 
certainty that smoothly approaches Heisenberg scaling as 
goes to zero ■ Due to their strong interactions, Bose- 
condensed atomic vapors are ideal systems for creating 
number-difference squeezed input states and using them 
as input into atom interferometers (lp| - [l3l l2p - [23| . For 
these systems, the GSS is an excellent approximation to 
the ground-state at T=0, with the squeezing controlled 
by adiabatic variation of the interaction-to-tunneling ra- 
tio 0. 

While strong interaction is essential for initial number- 
squeezed state preparation, it also limits the precision 
of the interferometer due to phase-diffusion during the 
phase acquisition time j25l - l33l ]. This process can be 
viewed as the shearing of the initial phase-space distri- 
bution due to the different mean-field shifts experienced 
at different points in the distribution. Since the mean- 
field shift is proportional to the population imbalance 
between the two condensates, phase-diffusion is propor- 
tional to the relative-number variance, A, during the hold 
time [H with characteristic decoherence time of 1/(UA), 
where U is the interaction strength. Number-squeezed 
states (A <C yN/2) are transformed to phase-squeezed 
states(A 3> s/N /2) by the first beam-splitter, thus pro- 
viding sub-shot-noise accuracy, at the cost of increased 
sensitivity to phase-diffusion. By contrast, states which 
are number-squeezed during the phase-acquisition pe- 
riod, are far more robust, but suffer from inherently large 
readout uncertainty. The interplay between readout un- 
certainty and robustness against phase-diffusion implies 
that the initial squeezing should be optimized to give the 
best possible precision [33[. 

It is conventional to describe the MZI as a device that 
measures 9, the path-length difference between the two 
arms of the interferometer. While this is the proper way 
to view an optical interferometer, trapped- atom interfcr- 
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ometers differ in that there is no fixed relation between 
time and distance. Thus we propose that the double- 
well condensate MZI be viewed as a device to measure 
the 'bias', or energy differential, e, between the two wells. 
Unlike 'flying particle' interferometers, the accumulated 
phase-shift 8 = eT, in a 'stationary particle' interferom- 
eter is not a fundamental measurable quantity, as the 
hold time T is a free parameter, which can be used to 
optimize the measurement of e. Here we use the free- 
dom of the initial number- difference variance, a, and the 
hold time, T, to optimize the signal-to-noise ratio of a 
Mach-Zehnder atom interferometer with a GSS input. 

The optimization is performed first using exact numer- 
ical results, and then using approximate analytic expres- 
sions, with the two approaches showing good agreement. 
As the optimal performance improves with decreasing in- 
teraction strength, we assume that the experimenter has 
reduced the collision strength, U , as much as possible, 
given the constraints of the experimental set-up. Our 
goal is then to proscribe the optimal squeezing and hold- 
time based on this minimum value of U. In the case 
where U « is obtained, e.g. via a Feshbach resonance, 
the experimental uncertainty, AU should be used in place 
of U. 

This approach connects previous work on optimizing 
a non-interacting interferometer [20j with the strong- 
interaction optimization of "useful squeezing" [33j, 
clearly showing a transition from optimal number- 
squeezing to optimal phase-squeezing, and mapping the 
transition region between the two regimes. 

The two-site Bosc-Hubbard model, initial state prepa- 
ration, and optimization function are presented in 
Sect. II. Numerical optimization results are described in 
Sect. Ill and compared to analytic predictions in Sect. IV. 
In Sect. V we compare the best accuracy obtained for 
Gaussian states with that attainable from a coherent in- 
put, with conclusions presented in Sect. VI. 



opposite 7r/2 readout rotation about J x (iii). The final 
population imbalance J/ is used to read the accumulated 
phase 8 = eT from which the bias, e, is readily obtained. 

Assuming that the beam-splitter and read-out rota- 
tions are instantaneous with respect to the characteristic 
phase-diffusion time, the MZI can be described by the 
propagator 

U M zi(0,u,j) = e _i * ■*» e -*-M 1 -(«/i)-?.) e - < * 3 % (2) 

where u = Uj/e. This propagator acts on a Gaussian 
squeezed state, of the form, 
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where a is the initial number-difference uncertainty, and 

M a = Em=-i ex P(~ m2 /( 2cr2 )) ~ %/2ttct- These states 
are an excellent approximation for the grou nd state of 
Hamiltonian ([1]) with U > and e ~ [24|], as well as 
to the dynamically squeezed states produced by single 
axis twisting of initially coherent states [H, [H, [TtJ , thus 
they can be readily generated with current experimental 
setups. 

Because the operators for phase-acquisition and phase- 
diffusion commute, the expectation value and variance J z 
at the output can be evaluated by applying the phase- 
diffusion operator to the state vector, while incorporating 
the bias via a rotation of the observables by angle 8 = eT. 
With the definitions 
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sin 6 J x — cos 8 J y , 
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II. MODEL AND INITIAL PREPARATION 



we see that (J z )out 
pressed as 



(v\UmziJz U mzi\<j) can be ex- 



We consider a Mach-Zehnder interferometer . realized 
via the two-site Bosc-Hubbard Hamiltonian [3J-l3a| 



H = -KJ X + eJ z + UJi 



(1) 



Here K, e, and U are coupling, bias, and interaction 
energies, where U > corresponds to repulsive inter- 
actions and vice versa. The bias e may be positive or 
negative, depending on the energy detuning between the 
two modes. The SU(2) generators J x = {a\a,2 + a\ai)/2, 
J y = {a\a,2 — a\a\)/(2i), and J z = (n\ — ri2)/2, are de- 
fined in terms of the boson on-site annihilation and cre- 
ation operators Sj, fit, with the conserved total particle 
number ni + n-i = N = 2j. The interferometer scheme 
(Fig. [TJa)) consists of a fast ir/2 beam-splitter rotation 
about J x (i), followed by relative-phase acquisition dur- 
ing a hold-time T due to the bias detuning e (ii) , and an 



(Jz)out = (<r,9\Jz,9\o;0) = sin6 {J x ) a< e, 



where 



(J, 



[i e {x, y, z}. 
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The uncertainty in the final measurement is then given 
by 



A ^z, out = sin2 & A Jx,a,8 + c °s 2 OA J, 



V,cr,t 



where 



AJ^,e = yJiJZUe-Wlf 



(8) 



(9) 



We note that there are no terms proportional to {J y ) a ,e, 
(J y Jx)cr,0, and (J x Jy) a ,0 in Eqs. (JB]) and © due to the 
symmetry of the state \<r,9). 
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FIG. 1: (Color online) Mach-Zehnder interferometry with 
squeezed input in the presence of collision-induced phase dif- 
fusion, (a) The interferometer sequence on the Bloch sphere, 
without interactions. An initial number squeezed state is (i) 
rotated by 7r/2 about J x by the first beam-splitter, (ii) ac- 
quires a phase difference during the hold time T, and (iii) 
counter-rotated about J x by the second beam-splitter to give 
a final Jz population imbalance readout, (b) Same with in- 
teractions. Phase diffusion results in the spreading of the 
squeezed states during the hold time, thus degrading the read- 
out accuracy, (c) Final number distribution and average pop- 
ulation imbalance (solid line) as a function of the acquired 
relative phase 9 for a non-interacting gas. (d) same in the 
presence of phase diffusion, demonstrating reduced fringe vis- 
ibility. 



It is useful to define the bias-measurement signal-to- 
noise ratio (SNR) as 



£ U 

From the error-propagation formula, 

n -1 



AO = 



d(J z ) 



out, 
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it is then straightforward to use ((SJ) and 



s(a,6) 



\0\ 



l J; T. 



AJ y 2 M + tan 2 0Aj2, 



(10) 

(11) 
to obtain 

(12) 



Because the state |<r, 9) depends only on the parameters 
{j, u, a, 9}, it follows that the optimal values, (co, $o) that 
give the maximum SNR, sq = s(ao,9o), as well as sq 
itself, are functions of j and u only. 



III. NUMERICAL RESULTS 

The goal of this paper is to find the optimal parame- 
ters, (a ,9 ), that give the maximum signal-to-noise ra- 
tio, s a = s(a ,9 ), for fixed U, j, and e; and then use 



this to determine the scaling of ,s with particle number, 
N = 2j, and interaction-to-bias ratio, u. We note that in 
the case s < 1, one should instead minimize the absolute 
uncertainty, Ae = A8/T. In this work, however, we will 
consider only the case s > 1, with the "minimum de- 
tectable bias", E m i n defined by s Q (j, u = Uj/e m in) = 1- 

The results of such optimization, using the numerical 
evaluation of Eq. are shown in Fig. [5] and Fig. [3] 
In Fig. [2] we plot the optimal squeezing <t , and the re- 
sulting maximized precision p Q = log 10 s (p directly 
corresponds to the number of significant figures of the 
read-out), as a function of the parameters u and j. In 
Fig. [3] we plot the optimal acquired phase 9 = eT a as a 
function of j and u. The fact that the acquired phase is 
optimized, rather than the hold time is discussed in de- 
tail in the next section. Symbols in Fig. [2] correspond to 
the analysis presented in Section IV. The optimal initial 
number variance increases with increasing interactions, 
whereas precision is degraded. At the limit of small u, 
we obtain that a scales as j 1 / 3 and s Q scales as ?' 2 / 3 
(dashed lines in Fig. [5]), in agreement with Ref. |20| . 
As the interactions increase, these power laws arc re- 
placed by a ~ \fj dependence of both quantities. Most 
significantly, a transition from optimal initial number- 
squeezing (a a < U j/2) to optimal initial phase-squeezing 
(<j > \J j/2) takes place as the interaction parameter 
u crosses unity. This transition results from the inter- 
play of projection noise minimization by initial number 
squeezing (i.e. hold-time phase-squeezing, resulting in a 
narrower 'phase-dial') and phase-diffusion control by ini- 
tial phase-squeezing (i.e. hold-time number-squeezing, 
rendering the state more robust against phase-diffusion). 
As the interactions increase, phase-diffusion determines 
the interferometer's precision and initial phase-squeezing 
is preferred. With current experimental set-ups, number- 
squeezed states are more readily obtained than phase- 
squeezed states [T3|. As the first-beam splitter rotates 
a number-squeezed state into a phase-squeezed state, 
initial phase-squeezing in the MZI picture is obtained 
in practice by preparing a number-squeezed state with 
(j — > and eliminating the first 'beam-splitter' J x - 
rotation, so as to have a number-squeezed state during 
phase- acquisition. 



IV. ANALYTIC OPTIMIZATION IN THE 
PRESENCE OF PHASE-DIFFUSION 

Noting that (J y ) a ,0 = identically and using Eq. ©, 
we may rewrite the signal-to-noise ratio of MZI output 
CCU as, 



s(a,6) = 



where 




(J 2 v 



(13) 
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FIG. 2: (Color online) Optimal values of the initial relative- 
number variance a and precision p. Left panels depict the 
dependence on the interaction parameter u at fixed atom 
number j = 64 (o), 128 (□), 256 (A), and 512 (o). Right 
panels show the dependence on j at fixed ln(it) = —2.5 (*), 
(+), and 2.5 (x). Solid lines are exact values whereas sym- 
bols correspond to the weak-interaction estimate of Eqs. (|32[1 
(filled) and the stron-interaction estimates of Eqs. (|36[) (un- 
filled). Dashed lines denote the j dependence for u — 0, 
whereas the dash-dotted line marks the width of a spin coher- 
ent state, separating initial number-squeezing below it from 
phase-sqeezing above. The transition from number- to phase- 
squeezing takes place at ln(-it) ~ 1 



The optimized signal-to-noise ratio is then obtained by 
minimizing Q(er, 9) with respect to <r for fixed 9 to obtain 
cr (9), and then maximizing s(a (9), 9) with respect to 9. 

In order to derive an approximate analytic expression 
for Q{<j, 9), we rely primarily on the approximation that 
a rotated Gaussian state is itself a Gaussian. Thus we 
make the anzatz, 



ln(u) 



FIG. 3: (Color online) Optimal acquired phase o for Gaus- 
sian squeezed states, as a function of the interaction parame- 
ter u and the particle number j. 



Replacing the sum by an integral with respect to x 
(m+l/2)/j gives, 
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u9 



for which A J, 



(15) 

= A. This variance is clearly a constant 
of motion during phase-acquisition (i.e. it is independent 
of u and 9). As the initial state |er) is a minimum un- 
certainty state with AJ Z = a, it follows that the state 
immediately after the first J^-rotation, |er, 9=0), is also 
a minimum uncertainty state, with AJ y = a. From the 
Hcisenbcrg uncertainty principle, it then follows that 



where we have made the approximation j + 1/2 « j, and 

_ j 2 

dropped terms proportional to e ^ . 

To proceed further, make use of conservation law, jj + 
J 2 + J 2 = J 2 w j 2 , which together with the relation 
Jl-Jl = \{J+ + J-), gives us 

^ODofi = f - A 2 + Re {(J 2 ) CT , e } , (19) 



and 



2(J*)„,e=f -A 2 -Re{(J 2 ) ff , e } 



(20) 



A = 



(Jx)a,6=0 

2a ' 



Thus in order to calculate Q(a,9), we need to compute 

(16) 

only (J x )a,e and (J?)^. Following the same procedure 



Noting that (J x )a,9=a = (<r\Jx\o) = H-7+k), with J± 
J x ±iJ y , we arrive at 



used to arrive at (fT5)l . we find 

1 2u 2 9 2 A 2 
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J exp 
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2er 2 V2^ 



E Vti- m )(j+ m + 1 ) 
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and 



j 2 exp 



2A 2 



8u 2 6 2 A 

~J 2 ~ 



A 2 

J 2 



16u 2 9 2 A 



(22) 



Combining these results, and expanding in terms of the 
three independent small parameters: 1/cr, cr/j, and the 
phase-diffusion winding- angle, 9 ( i — u9/cr; we arrive at 
Q(a, 9) = Q m (a, 9) + Qi(a, 0), where 



INI 



(a, 6) 



tan 



r 



32cr 4 ' 



(23) 



gives the effects of squeezing in the absence of interac- 
tions, and 



Ql(<7, 



1 + tan 2 



r 



(24) 



determines the effects of interactions. The extremum 
condition, do-Q(a,9) = 0, can be expressed as 



bx-1 



0. 



(25) 



where we have introduced the dimensionless parameters 
x = (vo/vni) 2 , b = (u/u c ) 2 , and a = ^/tan(0)/4, where 



• i /a 

(Tni = aj 1 , 
is the non- interacting solution (20| . and 

a 2 



(26) 



(27) 



is the critical collision parameter above which interac- 
tions predominantly determine the optimal performance. 

It is readily seen that in the non-interacting case (b — > 
0), we have x = 1, in which case the optimized signal-to- 
noise ratio is 



sni 



3 a 



(28) 



The optimal angle Oni is then found via (1 — 9de)a(9) = 
0, which has the solution 9^j = 1.14 radians, giv- 
ing a NI = 0.86J 1 / 3 , and s NI = 1.14j 2 / 3 . While the 
non-interacting MZI with GSS input can indeed exhibit 
"Hciscnbcrg scaling" at a = 1 and 9 = our optimiza- 
tion shows that it is more advantageous to instead use 
6 ~ 1 and a ~ j 1 / 3 ; a result that follows from the point- 
of-view that the property actually being measured is the 
bias, e, rather than the phase, 9. The Hcisenberg-limitcd 
measurement gives A9 = but with a signal-to- noise 
ratio of s ~ 1, smaller than s by a factor j 2 / 3 3> 1. 

Furthermore, this shows that in the non-interacting 
case, there is no minimum observable bias. A single 
measurement of any e, no matter how small, can yield 

f l°g'io J (meaning that 



a maximum precision of p 



increasing the atom number by x 30 results in one addi- 
tional decimal place of precision). In contrast, with no 
squeezing (er = y/J/2), the maximum obtainable preci- 
sion scales as \ log 10 j, and therefore requires an increase 
of x 100 atoms for each additional significant figure. Of 
course in practice, at long hold-times, phase-diffusion in- 
evitably degrades the performance, hence the optimized 
performance in the presence of interactions (u ^ 0) is ul- 
timately of more interest than the non-interacting case. 

In the weakly interacting regime b -c 1, Eq. ([23]) can 
be written as 



-l = bx. 



(29) 



Treating the r.h.s. as a perturbation gives the solution 

(30) 

resulting in the weak-interaction behavior, 
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The strong interaction regime !i > 1 was studied in 
Ref. [33|], using generic phase-space arguments for the 
dynamics rather than an explicit ansatz. The a 2 / j 2 and 
u 2 9 2 j a 2 terms on the r.h.s. of Eqs. (23) and (24) are re- 
spectively identical to the noise and phase-diffusion terms 
used in Eq. (20) of [33[ to minimize the "useful squeez- 
ing" , defined as ^/2]AJ v /\ (J x (a, 9)) \ [13, [H- However, 
assuming that the contribution of AJ X is negligible is 
only valid for initial preparations where neither A J* nor 
AJ l z significantly exceed the coherent-state variance y/j. 
This greatly reduces the range of a available for opti- 
mization. In particular when u = 0, ignoring the AJ X 
contributions to Q(a,9) fails to reproduce Eqs. (|2"6"|) and 
(|2"8"|) . In Sec. VI, we provide a detailed comparison of 
our results with those of Ref. [33| . 

In the strongly-interacting regime, b 3> 1, Eq. (f2"5j) can 
be expressed as 



x 2 -b= -. 
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Again, treating the r.h.s. as a perturbation gives 



resulting in 
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FIG. 4: (Color online) Optimal precision for a spin coherent 
state as a function of u at fixed j (left) and as a function of 
j at fixed u (right). Parameter values and notation are the 
same as in Fig. [2] with symbols corresponding to Eq. 1 )431 . 



To leading order, this gives the strongly-interacting re- 
sults, 



\puB~o 




The optimal acquired phase 9 a is obtained for weak 
and strong interactions respectively, by substitution of 
cr from Eq. ([32]) and Eq. ([36]) into d(Q/9 2 )/d0 = 0. 
Consequently, for strong interactions (u > u c ), 



29 2 a tan 9 Q (1 + tan 2 9 Q ) ps j/u, 
whereas for weak-interactions [u < u c ) we have, 

2 



sin 29 1 



2 u W/ 3 2 
tan i 



:9n 



(38) 



(39) 



When it = 0, Eq. ([39]) reduces to sin(20 o )/(20 o ) ps 1/3 , 
resulting in the appropriate non-interacting solution 9 a = 
8ni,- Eqs. ([3"5]) and (|5i?]) are in good agreement with the 
numerical optimization shown in Fig. [3] 

The weak-interaction Equations ([32]) . ([39 ]) and the 
strong- interaction Equations ([56"]) . ()38)) separated by the 
condition u = u c , constitute our main result. As shown in 
Fig. [2] they agree very well with the numerical optimiza- 
tion results. The weak-interaction power-laws a oc j 1 / 3 , 
s cx j 2 / 3 are continuously replaced by the scaling of 
both quantities as the interaction strength is increased. It 
is also evident from Eq. (]36p . that the optimal squeezing 
changes from initial number squeezing to initial phase- 
squeezing at u9 ps 1. 



V. COMPARISON WITH COHERENT INPUT 

It is instructive to compare the optimized SNR with 
that of an initial spin coherent state exp[i7T J y /2]\j, — j), 
approaching the a = W 'j jl Gaussian, in the presence of 
interactions. Numerical optimization results are shown 




(37) FIG. 5: (Color online) Same as Fig.[3]for spin coherent states. 



in Fig. [4] and Fig. [5] The phase diffusion of the coherent 
state is given as, 
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[1 + (cos2t) 2j '~ 2 - 2(cosr) 



^[l-(cos2r) 2 ^ 2 ] , (41) 

r [l-(cos2r) 2 -''- 2 ] 

4[l + (cos2r) 2 ^ 2 ] , (42) 



where r = (u/j)9. For u = the optimal phase is 
9 C = 7r/2, minimizing projection noise and resulting in 
the SNR s c = TVy/j/2. In the presence of interactions 9 C 
(Fig [5]) decreases to reduce the phase-diffusion time. The 
optimal SNR is given by, 



1/2 



s c = 0cV2j (l + 4?n w 
with the optimal relative-phase given by 



26» 3 tan 9 r ( 1 + tan 2 



.l/u 2 



(43) 



(44) 



The best SNR for a coherent preparation thus approaches 
s c PS yjj/{2u 2 ) for u8 c > 1. 

In Fig. [6] we compare the best precision obtained for 
spin coherent states (bold lines) with that of the op- 
timized Gaussian squeezed states. For weak interac- 
tions, the coherent states preform worse than the ini- 
tially number-squeezed states due to the larger projection 
noise, with the anticipated j 1 / 2 vs. j 2 / 3 respective signal- 
to-noise ratios. However, as the interactions increase, the 
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actions are sufficiently large. The transition from opti- 
mal number-squeezing to optimal phase-squeezing takes 
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FIG. 6: (Color online) Comparison of best precision for Gaus- 
sian squeezed states (normal lines) vs. spin coherent states 
(bold lines). Left panel shows the dependence on u at fixed 
j = 64 (dashed) and 512 (solid). Right panel depicts the de- 
pendence on j for fixed u = and as a function of j at fixed 
u = (dashed), ln(tt) = (dash-dotted), and ln(u) = 2.5 
(solid). 



coherent states are less affected by phase diffusion due to 
their smaller hold-time number variance. The number 
squeezing of the optimal Gaussian state decreases with u 
in order to slow down phase-diffusion, until at u = u c it 
coincides with the coherent state. Beyond this point, the 
optimal Gaussians are initially phase-squeezed, rotated 
to number-squeezed states during the phase acquisition, 
thus slowing down phase-diffusion with respect to the co- 
herent states. While retaining the same s cx j 1 ' 2 scaling 
at fixed u, the best Gaussian states offer a factor of ^fu 
improvement in the SNR over coherent states, as evident 
from comparison of Eq. ([36]) and Eq. ([43]) . 



place at u ~ u c cx j 



,-1/3 



The scaling of the best SNR 



with j changes from the sub-shotnoise interaction-free 
Eqs. ([26]), ([28]) [2(J through the weak-interaction Eq. ([32]) 
to the strong-interaction behavior of Eq. ([36l) . 

Before closing, it is instructive to compare in detail 
our results in the strong-interaction regime with those 
of Ref. [33|. This work optimizes the useful squeezing 
with respect to a for any given acquired phase 9. This 
approach amounts to neglecting the AJ X contributions 
and retaining only the a 2 /j 2 term in Eq. (|23[) and the 
u 2 9 2 /a 2 term in Eq. ([2~4"]) . The optimal a for any fixed 
value of 9 is a = \/Ju9 = jy/r, giving Q = 2u9/j = 2t 
and s — 9 j \J~Q~ = \j9jj (2u). These expressions are the 
same as those appearing in [33j and superficially seem 
to coincide with our Eq. (|3"6"]l . However, checking for 
self-consistency by substituting a Q into ([23]) and (|24|) . 



VI. CONCLUSIONS 



we obtain that (AJ X ) could be neglected with respect 
to (AJ y ) 2 only when j 1 / 2 < 2er < j 2 / 3 , because the 
A J x variance grows for both number- and phase-squeezed 
states, whereas the AJ y increases monotonically with 
hold-time number-squeezing. This greatly restricts the 
validity of the useful-squeezing optimization to values of 
9 for which a Q lies in the appropriate squeezing window. 
By contrast, our calculation extends to all values of 9 and 
a by including also the A J x contribution. 

In order to maximize the interferometer precision, it is 
imperative to devise schemes which will overcome or con- 
trol phase-diffusion. One such approach may be to sepa- 
rate phase-acquisition from phase-diffusion, by reversing 
the roles of e and K, so as to measure frequency-shifts 
of Rabi oscillations. This will allow for the use of robust 
number-squeezed state without loosing readout accuracy, 
at the expense of having to follow an essentially nonlinear 
oscillation. Future work will also seek to similarly opti- 
mize an SU(1, 1) interferometer, based on the stimulated 
dissociation of molecular BECs 1321] . 



Interaction-induced phase-diffusion is currently the 
most prominent obstacle in the way of sub-shotnoise 
atom interferometry [33|. By optimizing the SNR of a 
Mach-Zehnder atom-interferometer in the presence of in- 
teractions, we confirmed and quantified the notion that 
whereas phase-squeezing during the phase-acquisition 
time is required to go below the standard quantum 
limit, the robustness of number-squeezed states to phase- 
diffusion makes them the preferred choice when inter- 
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